Far from equilibrium quantum magnetism with ultracold polar molecules 
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Recent theory has indicated how to emulate tunable models of quantum magnetism with ultracold 
polar molecules. Here we show that present molecule optical lattice experiments can accomplish 
three crucial goals for quantum emulation, despite currently being well below unit filling and not 
quantum degenerate. The first is to verify and benchmark the models proposed to describe these 
systems. The second is to prepare correlated and possibly useful states in well- understood regimes. 
The third is to explore many-body physics inaccessible to existing theoretical techniques. Our 
proposal relies on a non-equilibrium protocol that can be viewed either as Ramsey spectroscopy or 
an interaction quench. It uses only routine experimental tools available in any ultracold molecule 
experiment. 
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Excitement about the recent achievement of near- 
degenerate ultracold polar molecules [THE] in optical lat- 
tices [6] stems from their strong dipolar interactions 
and rich internal structure, including rotational, vibra- 
tional, and hyperfine states. These features may be 
applied to tests of fundamental constants [7], quan- 
tum information [8 , ultracold chemistry [9], and quan- 
tum emulation of condensed matter models [T0HT2] . In 
this paper our focus is on molecules as emulators of 
quantum magnetism [13-24 , specifically as proposed in 
Refs. [25j|26]. Models of quantum magnetism have some 
of the simplest many-body Hamiltonians, yet describe 
numerous materials [27U29] and display condensed mat- 
ter phases ranging from fundamental to exotic: anti- 
ferromagnets, valence bond solids, symmetry protected 
topological phases, and spin liquids. Emulating quan- 
tum magnetism with molecules is appealing because, like 
cold atoms, the systems are clean and the microscop- 
ies well understood. Advantages over cold atom em- 
ulations of quantum magnetism [30 include orders of 
magnitude larger energy scales and more tunable Hamil- 
tonians [25j [26]. These prior studies have focused on 
spin ground states of unit filling insulators. In con- 
trast, we propose a simple dynamic procedure applicable 
to present experiments, which are ultracold, but non- 
degenerate and low density. We show that interesting 
many body quantum magnetism can be studied immedi- 
ately. 

Specifically we show how experiments may use this 
dynamics to achieve major goals for emulating quan- 
tum magnetism, and we outline these goals to motivate 
our calculations. First, although interesting models of 
quantum magnetism are predicted to describe ultracold 
molecules under appropriate circumstances, this has yet 
to be experimentally demonstrated. The proposed dy- 
namic protocol allows such a demonstration as well as 
benchmarking of the emulator's accuracy. Second, one 
wishes to prepare interesting correlated — and possibly 
useful — states. This protocol can generate such states 
in well-understood regimes. Finally, one wants to ex- 




FIG. 1: (Color online) Dynamic protocol viewed as Ramsey 
spectroscopy, two microwave pulses of area and 0' sepa- 
rated by time t. Inset: alternatively viewed as an interaction 
quench, a sudden t = turning off of an infinitely strong field 
he along 0. 



plore behavior in these models in regimes inaccessible to 
present theoretical tools. This is the generic case for the 
proposed dynamics. We emphasize that all of these goals 
are achievable under existing experimental conditions [6] , 
despite present experiments being non-quantum degen- 
erate and at low density. Furthermore, they require only 
routinely used measurement and preparation tools [3T] . 

Background. — Refs. [25, 26 show how molecule rota- 
tional states can serve as effective spins, and that dipolar 
interactions provide an effective spin-spin interaction. In 
the simplest case, one populates two rotational levels in a 
dc electric field E [67 j and works in a deep lattice to allow 
no tunneling. In this limit, a spin- 1/2 dipolar quantum 
XXZ model describes the molecules l68l: 
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The sum runs over all molecules, Sf and are the spin- 
1/2 operators satisfying [S?,Sf] = ±5^, and Vdd(hj) = 
(1 — 3 cos 2 ®ij)/\vi — Tj\ 3 with the z'th molecule's po- 
sition in lattice units and the angle between E and 
— Yj. For simplicity and concreteness we assume a 
dimension d < 2 system with E perpendicular to it, so 
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Vdd(h j) = Vl r i — r j| 3 ? but our ideas apply in arbitrary 
geometries. One may tune J±/J z by changing E and 
the choice of rotational state. We denote by |0), |1), 
and 1 2) the three lowest energy rotational eigenstates in 
the applied E-field with zero angular momentum along 
the quantization axis. Choosing |0) and |1) to make the 
spin- 1/2, one can tune oo > J±/J z > 0.35 using (readily 
achievable) E- fields from to 16 kV/cm. Choosing |0) 
and 1 2) for the spin-1/2, one can tune < J±/J z < 0.1 
for similar E- fields. A characteristic scale for these cou- 
plings is 400 Hz in KRb and 40 kHz in LiCs [32], com- 
pared to <10Hz in cold atoms using superexchange [33] . 
KRb molecules recently have been loaded in a deep three- 
dimensional lattice with 25 second lifetimes [6], allowing 
dynamics lasting thousands of JJ 1 and J z x . 

One important aspect of the ongoing experiments is 
that the filling / is much less than one molecule per site. 
The JILA experiments estimate / ~ 0.1. As a result, 
although the molecules' positions are static throughout 
one shot, they fluctuate shot-to-shot. Thus, rather than 
forming a regular lattice, the spins' locations have signif- 
icant disorder. Our calculations show that the dynamic 
protocol's utility persists with disorder. 

We use a simple disorder model that likely describes 
current experiments. We assume that each site is occu- 
pied with a probability p that is independent of other 
sites [69 . If the molecules are fermions (e.g., KRb [6]) 
then for current temperatures, which occupy only the 
lowest band, no sites can be doubly occupied and p = f. 
This also applies to bosons with a strong on-site density- 
density interaction (e.g., RbCs [34]). The trap causes / 
to vary spatially. Although we show results only for the 
homogeneous system, we have taken the trap into ac- 
count and found that our conclusions remain valid [35] . 

Remarkably, close relatives of such seemingly unusual 
models exist in the literature, for example the Blume- 
Emery- Griffiths model [36]. These mainly focus on the 
classical equilibrium limit, J± = 0. They were intro- 
duced to understand materials [37J [38], 3 He- 4 He mix- 
tures [36] , Griffiths phases [39] , glassy dynamics [40j [41] , 
and the interplay of the underlying lattice's statistical 
mechanics (site-dilution percolation) with that of the 
magnetism living on that lattice [42-44 . Rather than 
studying unique disorder-induced behavior, we focus on 
showing that / = 1 behavior survives disorder. 

Only through "disorder" does temperature enter, be- 
cause the deep lattice freezes out the motion. In particu- 
lar, one must distinguish motional temperature from spin 
temperature. The former may be large but is entirely 
captured by the disorder, while the latter is ill-defined 
since we consider non-equilibrium spin states. However, 
experimental microwave manipulation can produce es- 
sentially zero entropy spin states. While one could worry 
that disorder washes out the behavior, we will show that 
strong correlations, entanglement, and interesting many 
body physics survive large amounts of disorder. 



Dynamic protocol — Our dynamic procedure may be 
alternatively viewed as Ramsey spectroscopy or an in- 
teraction quench (Fig. [I]). Ref. [48] studied closely re- 
lated Rabi spectroscopy. In Ramsey spectroscopy, a well 
established tool in atomic physics, one begins with all 
molecules in the rotational ground state and applies two 
strong, resonant microwave pulses separated by time t. 
The first pulse initializes the spin states along 0, specif- 
ically to cos(<9/2)e^/ 2 ||) + sin(<9/2)e-^/ 2 |t), for an an- 
gle set by the pulse area, with high fidelity (> 99%). 
We take ip = with no loss of generality. The sec- 
ond pulse rotates a desired spin component, chosen by 
the pulse area and phase, to the z axis. In this way 
one can measure any desired collective spin component 
(n • S), where n is a unit vector and S a = Sf with 
a G {x,y,z}. One can also obtain higher moment cor- 
relations, e.g. ((n • S) 2 ), from the measurement record. 
Between these pulses the spins evolve for a time t under 
the Hamiltonian in Eq. ([I]) . We note that molecule exper- 
iments have recently begun using this protocol [31] and 
Ref. [45 applied it to long-range Ising models in recent 
Penning trap experiments with ~ 300 ions. 

If one imagines adding a transverse field term hS • he to 
the Hamiltonian in Eq. ([!]), with he a unit vector point- 
ing 6 from the — z axis (see Fig. [TJ, the Ramsey protocol 
corresponds to a quench from h = oo to h = 0. One 
may therefore be able to explore, for example, Kibble- 
Zurek physics (e.g., entropy production, topological de- 
fects) 051133. ' 

Theoretical methods. — We calculate dynamics in four 
limits: (1) short times, { Jj_, J z }t <C 1, (2) Ising, J± = 0, 
(3) near-Heisenberg [SU(2)j, \J Z - J±\ <C J z , and (4) 
one dimension for arbitrary J±/J z . The first three 
limits' results are analytic and valid in any dimension, 
while the last is from numerically exact adaptive time- 
dependent density matrix renormalization group (adap- 
tive t-DMRG) [49U52] . Details of the calculations will 
be presented elsewhere [35]. In all cases S z is conserved, 
with (S z (t)) = -(//2)cos0. 

Short time limit, {J_\_,J z }t <C 1. For short times, 
(0(t)) = (0)-it([0,H])-^([[0,H},H})+0(t 3 ) for 
an operator O. We calculate the commutators and time 
dependence of (S a (t)} to leading non-zero order, and 
(S a (t)S^(t)) to linear order. We find 

(Sf) = ^sin^|l-^[S 2 + /Tcos 2 ^]}+0(r 4 ), 
(S?) = -(ASx/Sj^flJ + OCr 3 ), (2) 

where r = (J z - J±)t, E m = 5^ Kid(M + j), and 
T = S 2 — S2. Note that for these homogeneous systems, 
these observables are independent of i. Similarly, defining 
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FIG. 2: (Color online) "Phase diagram" illustrating 
crossovers of dynamics versus filling /, J±/J z , and using 
t-DMRG on chains. In each region, a plot shows dynamics 
for = 0.l7r,0.57r (top blue, bottom purple), labeled with 
a qualitative description of the behavior. Shaded regions in 
/ = 0.2 plots indicate one standard deviation errors from dis- 
order averaging. "KRb: 14 ms" indicates the time for KRb 
in a one-dimensional 532 nm chain at a 5 kV/cm dc electric 
field giving J±/ J z — 2 for the two lowest energy, zero angular 
momentum projection rotational states. Dynamics in higher 
dimensions is faster due to having more neighbors. 
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(S?S]) : we find 

r/ 3 sin(2<9)sin<9 
16 

r/ 3 rsin(2<9)cos6L Tr ,. . 3n 
Si + V dd (i,j) sm 3 



[F dd (z,j)-Si]+0(r 2 ) (3) 
+ 0(r 2 ) 



for i j. To linear order, Cf^ and Cff are constant. 
For i = j, the Pauli algebra reduces (Sf Sf) to (Sf). 
One can compute E m rapidly for arbitrary lattices and 
analytically in special cases (e.g. one dimension). 

Ising limit, J± = 0. We extend the Emch-Radin solu- 
tion [53H55] for Ising dynamics to arbitrary inter-spin 
coupling strengths, and to include disorder [70 . We find 



(sm) = /^Re 



nE^ft) e ^ Jsl/,ld( '' ,+,)ft 

j^i Pi 



,(4) 



where the sum runs over p = (unoccupied site) and 
p= ±1 (S z = ±1/2), and 



g(p) 



i - / if p = o 

fsin 2 (<9/2) if p = 1 
[/cos 2 (0/2) if p = -1 



(5) 



The expectation (Sf) takes the imaginary (rather than 
real) part of the square-bracketed expression in Eq. Q. 
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FIG. 3: (Color online) Short time behavior as functions of 
6 and / in a one dimensional dipolar chain. Left: coefficient 
A/f defined by (Sf(t)) = (Sf (0)) - At 2 + 0(t 4 ) for fillings 
/ = 0, 0.2, . . . , 1.0, bottom to top, with r = (J z -J±)t. Right: 
B, C (top, larger and smaller curves, respectively), and D 
(bottom), defined as (S?) = -Br, (SfS?) = C(V dd (iJ) - 
Si)r (Si defined in text), and (Sf Sj) = Dr, respectively (to 
linear order). The coefficient D is shown for i — j = l,2,3,oo 
(top to bottom). Results are from Eqs. |2]) and (|3j. 



Similarly one can obtain correlations [35 , 56]. The prod- 
uct in Eq. Q is readily evaluated numerically by trun- 
cating the interaction range, even for a truncation in- 
cluding thousands of sites. In special limits (Sf) sim- 
plifies: e.g., for = tt/2 and / = 1, (Sf(t)) = 
(l/2)n^cos(J,y dd (i,j)t/2). 

Near-Heisenberg limit, J z « J±. — Here a finite size 
gap A oc J±/N 2 for N particles to excitations out of 
the Dicke manifold prevents states initially in the man- 
ifold from leaving it [53 [58] . Since matrix elements to 
leave it are ex J z — J±, projecting to the Dicke manifold 
is accurate when \J Z — J±\ « A. In this limit the effec- 
tive Hamiltonian is the collective spin N/2 model [57] [58] 
#eff = X(S Z ) 2 with X = iffw^ Ei^j V dd (i, i + j). Dy- 
namics are straightforwardly calculated for any disorder 
configuration, since there are only N + 1 states in the 
Dicke manifold. Here we restrict to / = 1 for simplicity. 
For example, one finds 



(Sf) = ^Re 



(cos(xt) — i cos 6 sin(xr)) 



N-l 



(6) 



Again (Sf) is the corresponding imaginary part. Unlike 
the other approximations, this is valid only for finite N. 

One dimension. — We use adaptive t-DMRG [49H52] to 
calculate dynamics of one dimensional chains. We treat 
20 site chains and find finite size effects to be fairly small. 
We discretize time in steps of 0.05J" 1 , and find a dis- 
carded weight of < 10 -9 for times < 10 J" 1 , adaptively 
keeping m = 50-500 reduced density matrix states. Al- 
together, we expect errors dominated by the disorder av- 
erage, which is taken over 100 random configurations. 

Results: global perspective. — Fig. [2] overviews dynam- 
ics, from the calculations above, as a function of J±/J z , 
/, and 0. Experimentally, these are controlled by elec- 
tric field [25j|26], temperature/density, and first Ramsey 
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FIG. 4: (Color online) Correlation dynamics and comparison 
with ground state: (S* Sj +i ) - (Sj) (Sj +i ) (averaged over j) 
for J± = 2J Z with ^ = 7r/2 as a function of time for / = 0.4 
(left) and / = 1 (right), compared to the ground state (upper 
right bar). Other /, 9, and J± are similar. 

pulse area, respectively. Fig. [2] shows dimension d = 1 
results, but our analytic expressions show that the d = 1 
results are representative of d > 1. Dynamics in d > 1 
have more neighbors and thus is faster. 

Consider / = 1 and = tt/2. For J ± /J z = 0, (Sf) os- 
cillates with period 2tt/J z from the nearest neighbor in- 
teraction, superposed with slower oscillations from longer 
range interactions. The first- few- neighbor interactions 
account for the dynamics to times t ~ 10 J" 1 . For 
J± = the frequencies form a discrete set. Increasing 
J± gives a continuum of frequencies, damping the oscil- 
lations. Approaching J_\_ = J z , the dynamics slows down, 
since at J± = J z the initial state is an eigenstate of the 
Hamiltonian. As J±/J z increases further, the dynamics 
is damped with characteristic timescale (J± — J z ) _1 . 

For / < 1, the behavior crosses over to that of in- 
dependent clusters, eventually with only two particles. 
The largest frequency is roughly half that for / = 1, 
since there is a single neighbor instead of two. Thus, the 
dynamics remains roughly as fast as for / = 1, but the 
dynamics' magnitude at times ~ 1S smaller 

since there are fewer molecules and only a fraction of 
them are close enough to interact. At any / the overall 
timescales are roughly independent of 6 but the damping 
vanishes as — » 0. 

Achieving goals of emulating quantum magnetism. — 
Fig. [3] shows the characteristic dependence of the XXZ 
model's short time dynamics on 0, /, and {J Z1 Jj_}. This 
can be used to experimentally verify the emulation of the 
XXZ model and benchmark its accuracy. 

To achieve the second goal of generating interesting 
well-understood states, both the near-Heisenberg and 
Ising limits are useful. For / w 1 and = tt/2 near the 
Heisenberg point, the state at t = tt/(2x) is \GHZ) = 

(l/\/2)(H <-) + |-> ►» for some EMQ]. 

This is a cat state, specifically the GHZ or NOON state, 
which is useful for metrology [61] . Ising dynamics offer 
other interesting states. For nearest neighbor interac- 
tions and = tt/2, the state at t = tt/(2J z ) is a clus- 



ter state, which suffices for universal measurement based 
quantum computing [62] . The presence of long range in- 
teractions perturbs the cluster state, and an interesting 
question is how this affects its utility. Decoherence can 
also limit the creation of entangled states. 

A generic implementation of the proposed dynamics in 
d > 1 achieves the third goal, emulating quantum mag- 
netism in theoretically intractable regimes. Away from 
the short time, Ising, and Heisenberg limits, no solution 
is known in d > 1. As Fig. [4] shows, in d = 1 strong 
correlations develop, suggesting the difficulty of d > 1 
calculations. The long time / = 1 correlations are even 
larger at large distance than in the ground state. Inter- 
estingly, the dynamics shows a light- cone-like spreading 
to an apparent steady state. 

Experimental outlook. — Though our discussion focused 
on molecules, we point out that the dynamics stud- 
ied here can have direct application in other physi- 
cal systems, including condensed matter [41 , trapped 
ions [45] [63] , and optical lattice clocks [64] [65] . 

We close by noting technical details for molecule exper- 
iments. Rotational states' polarizabilities differ [25] [66], 
so the optical trap induces a spatially varying field 
^ZihiSf. Also, Eq. ignores density-density n^nj and 
density-spin n^J interactions [25] [26]. For / < 1, the 
latter gives a spatially varying magnetic field that de- 
pends on molecules' random positions. Spin-echo pulses 
common in Ramsey experiments remove both effects. 

Summary. — We have shown that Ramsey spectroscopy 
enables ongoing ultracold polar molecule experiments 
to accomplish three goals for emulating quantum mag- 
netism: (1) benchmarking the emulation's accuracy (us- 
ing short time dynamics), (2) generating strongly cor- 
related and entangled states in well-understood limits 
(Ising, near-Heisenberg, one dimension), and (3) explor- 
ing strongly correlated dynamics in regimes inaccessible 
to theory (generic case in dimensions d > 1). 

Finally, we mention that in addition to the XXZ Hamil- 
tonian explored in this paper, our dynamic protocol 
should be useful for verifying emulation of more com- 
plicated spin models that may be realized with ultracold 
molecules, as in Refs. [25] [26] and beyond. 
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